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FILTRACION A PRESION DE PASTES CERAMICAS. IV - 
TRATAM1ENTO DE LOS DATOS EXPERIMENT ALES. 
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RESUMEN / 

Dado el elevado numero de tests que es neccsarid rcalizar para l^letcrminacion del comportamien- 
to de las distintas barbotinas en el proccso de filtracion a presioiyrc considcro intcresante la utilizacion 
de nuevas tecnicas de tratamiento dc dato que con un numcrc/rcducido de datos experimentales per- 
mitieran la realizacion de un ajuste riguroso del parametro nyoc la ecuacion representativa del proccso, 
directamcnte relacionado con la rcsistcncia cspeci'fica dc laiorta. 

En el presente trabajo se expone la utilizacion de un/metodo de tratamiento de datos basado en la 
aplicacion de algoritmo propuesto por Kalman, al prooCsd de filtracion a presion constante. asf como 
sus ventajas. Por ultimo se realiza una comparacion con el metodo de los minimos cuadrados. 

SUMMARY 

As the study of the behaviour of different kinds of barbotines in pressure filtration requires a con- 
siderable number of tests, the convenience of using new developed data processing methods was stu- 
died, which, based on a small quantity of experimental data, would allow the stringent adjustment of 
the parameter of the equation representing the process, which is in direct relationship to the specific 
strength of the lining. 

In the present paper the use of a data processing method based on the algorithm proposed by Kal- 
man and its application to the filtration process at constant pressure arc described as well as the advan- 
tages of this method. Finally this tecnique is compared with the method of the smallest squares. 

RESUME. f) 

Etant donnee la prande quantite de tests necossaires pour determiner le comportement des diffe- 
rentes barbotines pendant del proces de filtration a pression. on considere interessant le recours a des 
tecniques nouvelles d’utilisation de donnees quf permettraient avec une quantite reduite de donnees 
experimentales, la rdlisation d’un ajustege ripo/reux du parametre (a) de i’equation representative du 
proces, directement relationne avec la resistant cspecifique dc la Tourte. 

Dans ce travail on explique 1’utilisation dame methode de traitement de donne’es qui a pour base 
l’aplication de 1’algorithme propose par Kalmin, au proces de filtration a pression constante. aussi bien 
que ses avantages. En fin on fait une comparison avee la methode des minimes carrcs. 

ZUSAMMENFASSUNG. * » 

In Anbetracht der grossen Anzalil von Versjuchen. die zur Bestimmung des Verhaltens versciiiede- 
nartiger Schlichker wahrend des Druckfiher p/ozesses erfordcrlich sind. wurde der Einsatz neuartiger 
DV-Verfahren als lohnend betrachtet. die mi/einer geringen Menge expcrimenteller Daten die exakto 
Einstcllung des Parameters' der den Prozess durstellenden Gieichung gestatien. der zu der spczifischen 
Festigkeit des Kitchens in Bezienhung steiit. / 

In der vorliegenden Arbeit werden eineiuuf der Benutzung des von Kalman vorgesehlagenen Aigo- 
rithmus beruhende Methode und ihre Anwendung auf den Filterprozess bei konstantem Druck bescit- 
rieben; die Vorteile dieses Verfahrens wejuen erlautert. Abschliesscnd erfolgt ein Vcrgleich dessdben 
mit der Methode der Klcinsten Quadrate. 


1. INTRODUCCION. 

Para determiiiar ios parametros caracterizadores de una 
barboiina en el proccso de filtracion a presion y predecir su 
comportamiento a escala industrial, es necesario realizar un 
elevado numero de tests de laboratorio. por lo que resulta 
ntuy interesante la utilizacion de tccnicas de tratamiento de 
datos lo mas eficientes posiblc que proporcionen un aliorro 
del numero de pruebas a rcalizar usf como un atuilisis esta- 
distico de los resultados obtentdos. 

En el caso que nos ocupa la mayorta de los trabajos reu- 
lizados para mejorar la planificaeion de los experimentos v 
el tratamiento de los datos se h an orientado liacia el ctilciiln 
de los valores de la rcsistcncia cspeci'fica (o ctialqnier olro 
parametros relacionado con ellu) para predecir los resuiia- 
dos a escala industrial, f 1 , 2). 


En este trabajo se propone un procedimiento Je trata- 
miento de datos experimentales, ya descrito en la bibliogra- 
tVa, para la obtencion dc dicito parametros y de su fiabili- 
dad. en funcion de los datos utilizados. 

1.1. ECUACION TEORICA. 

Como se indico en la primera parte de este trabajo i 1 .). la 
ecuacion teorica representativa de la filtracion a presion es: 

t ppx.ee p R m 

— = V + ( !. ) 

V ( ! -M.x)2k.AJ\A 2 k.&l'.A 

qnc puede escribirse en la forma simplificada 
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‘ - 1 = a V + b (2) 

V 

Dado que en la mayoria dc los casos h resistencia ofreci- 
da por el medio filtrante es despreciabie frcnte a la resisten- 
cia de la torta ( 1 ,2 3 »), la evaluation del coefidente b rcsul- 
ta innecesaria. Por el contrario, la determination del coefi- 
ciente.a es de gran importancia porserelque caracteriza la 
filtrabilidad de una barbotina dada, para unas condiciones 
determinadas. 

2. ALGORUMO DE KALMAN. 

El procedimiento propuesto es el resultado de la aplica- 
cion del algoritmo de Kalman (4 y 5) al proceso de filtra- 
tion a presion constante. 

Este algoritmo permite ir obteniendo unos valores “me- 
jorados” de las variables y parametros, cada vez que se in- 
troduce un nuevo punto experimental, uttlizando un mode- 
lo numerico preestablecido. Para ello se <emplea la expre- 
sion: 



Estimation mejorada 
para el tiempo i. 


Como puede observarse, en este caso el modelo consta 
de dos variables y un solo parametro. Por lo tanto, unica- 
mente se podra obtener informacion acerca del valor de la 
pendiente de la recta (parametro a) ya que el calculo de ca- 
da punto se realiza en base al anterior en el tiempo, er. el 
cual queda implicito el valor de la ordenada en el origen. 

Hay que tener en cuenta que dado que tanto cl volumen 
como el parametro a no son valores exactos sino que vienen 
afectados de error se representaran mediante su media y su 
covarianza: 

Vj 'V (Vj , -2 vp 

a j ^ (*j > I a j) 

Asimismo hay que considerar el error asociado al modelo 
numerico (ecuacion (4) ), expresado por wj 'V(0.Q);es de- 
cir, se considera con valor medio cero y covarianza Q cons- 
tante. 


Correction debida a 
las nuevas medidas 
realizadas en el tiem- 
po i. 


h i/j 

Prediccion basada en 
la informacion dis- 
ponible hasta el 
tiempo j (i >j). 


donde 

h: vector formado por las distintas variables que intervie- 

nen. 

fj; vector diferencia entre los valores predichos y los me- 
didos para un tiempo i. 

Gj: matriz de ganancia de Kalman, que determina la in- 

fluence relativa de las componentes del vector 

Kalman utiliza como parametro estadi'stico caracteriza- 
do de la exactitud del me'todo el valor medio del cuadrado 
de los errores (diferencia entre el valor “mejorado” y el va- 
lor exacto), y trata de minimizarlo. 

A continuation se expone la aplicacion de este metodo 
para el tratamiento de los datos experimentales obtenidos 
en el proceso de filtration a presion constante y que fueron 
presentados en trabajos anteriores. 

2.1. MODELO MATEMATICO. 

La ecuacion (2) puede reescribirse, definiendo una nueva 
variable 1 =t/v, en la forma: 

1 = a V 4-b (4) 

El error de modelacion introducido al tontar como varia- 
ble el cociente t/v , la no constancia del parametro a, a lo 
largo de todo el proceso de filtration a presion (2.6), as! co- 
mo las distintas suposiciones hechas para tomar la ecuacion 
(4), como representativa del proceso, se tienen en cuenta en . 
este metodo, mediante la introduction, en el modelo pre- 
dictivo, de un error global asociado al modelo. 

Tomando V como variable dependiente y poniendola en 
funcion de su valor en el tiempo anterior, Vj queda: 

Al 

V = V, 4 — . (5) 


2.2. APLICACION DEL METODO. 

La aplicacion de este me'todo incluye dos pasos. 

En el primero es estima un nuevo valor del volumen y su 
covarianza para el tiempo i, a partir de la informacion obte- 
nida hasta el tiempo anterior. Este paso se denomina mode- 
lo de prediccion. 

A continuacion, partiendo del valor predicho para V y su 
covarianza y del valor experimental medido para ese mismo 
tiempo se recalculan a, V y sus covarianzas. Esta parte es 
propiamente la aplicacion del algoritmo de Kalman (ecua- 
cion 3). 

2.2.1. Modelo de prediccion. 

Para esta parte del calculo los valores del parametro a y 
su covarianza son los obtenidos en el intervalo anterior. Sin 
embargo, la incertidumbre del volumen, medida por su co- 
varianza, debe modificarse puesto qu'' es funcion no solo de 
la incertidumbre de V y a en el tiempo anterior, sino tam- 
bien del error asociado al propio modelo numerico. 

En el calculo de la nueva covarianza interviene la deno- 
minada matriz de sensibilidad que esta formada por las deri- 
vadas parciales de cada una de las variables respecto de cada 
uno de los parametros. En el caso presente, dicha matriz 
queda reducida a: 

3 

D= (Vj/:) (6) 

3a 

El valor de esa matrix, da idea dc la rapidez y bondad de! 
ajustc que sc esta reali/.ando. 

El calculo del valor predicho del volumen se realiza me- 
diante la aplicacion de la ecuacion 5. 
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i 2.2.2^ Aplicacion del algoritmo. 

• No se ha considerado convenierite exponer en este traba- 
jo el desarrollo promenorizado de las ecuaciones que cOnsti- 
tuyen el metodo propuesto, dado que estan suficientemente 
exph'citas en los trabajos originales. 

Para su aplicacion se ha desarrollado un programa BA- 
SIC, que es en este caso de una gran simplicidad y cuyo lis- 
tado puede remitirse a toda persona interesada en el. 

3. APLICACION. 

Como ejemplo ilustrativo del mdtodo propuesto se mues- 
tra el ajuste de los datos obtenidos en uno de los experi- 
mentos de filtracion a presion a escala de laboratorio. 

La barbotina utilizada fue porcelana y se trabajo a una 
presion de 8 Kg/cm 2 . 

Para la realization del ajuste se han hecho las considera- 
ciones siguientes: 

- el valor inicial de las covarianzas del volumen y del 
parametro se tomaron de tal forma que permitieran abarcar 
los intervalos de variation previsibles, para ambos. Los valo- 
res asignados fueron: 


|v = 100.00 
2 a = 0.001 

— el valor inicial de la covarianza cruzada, 2 v.a, se ha 
supuesto nulo, dado que no se conoce a priori el signo de la 
correlation. 

- el error de medicion del volumen se ha tornado igual a 
2 cm 3 intentando considerar no solo el error del aparato si- 
no tambie'n los errores sistematicos. 

A fin de comprobar la aplicabilidad del metodo propues- 
to se han realizado dos ajustes, para el mismo experimento, 
tomando en cada caso un niimero diferente de datos. 

En teori'a , el niimero y calidad de los datos experimen- 
tales necesarios para evaluar los pardmetros que definen un 
proceso son funcion de los valores de la matriz de sensibili- 
dad, que actiia a trave's de las covarianzas cruzadas y de la 
matriz de ganancia de los parametros. Un valor grande de 
esta matriz, provoca un ajuste rapido del parametro asf co- 
mo una diminution en la incertidumbre. 

En el caso en estudio la matriz de sensibilidad presenta 
un valor bastante elevado, por lo que en principio el nii- 
mero de observacior.es a realizar no ha de ser muy grande. 

El niimero de datos utilizados para los ajustes fue de 20 
y 8 respectivamente, tomando en ambos casos valores a lo 
largo de todo el proceso de filtracion. Los datos correspon- 
dientes a los tres primeros minutos de la operation no fue- 
ron considerados, ya que es necesario un cierto tiempo para 
establecer la presion de trabajo. 



Figura 1. 



Figuia 2. 
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Figura 3. 


En las figuras 1 y 2 se ha representado el valor del para- 
metro a y el de dos veces su desviacion tfpica. frente a 1 pa- 
ra los dos ajustes realizados. Puede observarse en ambos ca- 
sos una buena convergcneia. asf como un grado de incerti- 
dumbre suficientemente pequcno (del orden de centc'simas 
de la magnitud de a). 



En las figuras 3 y 4, se ha representado la diferencia en- 0 ° • 

tre el valor del volumen observado y el calculado asf como 

•el doble de su desviacion tfpica en funcion de 1. En ambos . > . • ' 

casos se puede observer que el error permanece en todo 

momcnto inenor que cl doble de la desviacion tfpica, valor Figura 4. 
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qu^ estSdisticamente represents el 90 /o de probabilidad 
de obtener el valor exacto. Dado que la disminucion cn cl 
grado de incertidumbre obtenido al realizar cl ajuste con 20 
valores en lugar de con 8 es muy pequena. puede conside- 
rate que es suficiente la realizacion de 8 6 10 observaciones. 
durante el proceso. 

Tambien hay que seiialar la gran simplicidad del metodo 
desarrollado, lo cual da lugar a un ahorro en el proceso de 
tratamiento de los datos. 

Por ultimo y a fin de comparar el me'todo propuesto 
con otros, se ha realizado un ajuste por el metodo de los 
minimos cuadrados. Los valores de a obtenidos mediante 
ambos metodos se recogen en la tabla I. pudiendose obser- 
var que las diferencias son muy pequenas. Pero hay que des- 
tacar que el mdtodo de los minimos cuadrados no suminis- 
tra ninguna informacion directa acerca de la incertidumbre 
del valor del parametro ni considera los errores de medida 
de las variables y los errores asociados al modelo matemati- 
co representativo del proceso. 

TABLA I 


n° puntos 


Metodo 

utilizados 

a 

Kalman 

8 

6.372.10" 3 

M.M.C. 

8 

6,385.10' 3 

Kalman 

20 

6.377.10' 3 

M-M.C. 

20 

6/379. 10" 3 


Tambien hay que tener en cuenta que un estudio estadis- 
tico muestra la jiecesidad de un mmimo de 30 valores a fm 
de poder aplicar la teoria de peuqenas muestras, y poder 
considerar la ecuacion obtenida mediante el me'todo de los 
minimos cuadrados como representativa del proceso estu- 
diado. 

De todo ello se deduce la importancia del algoritmo de 
Kalman en el tratamiento de los datos del proceso de filtra- 
tion. 

CONCLUSIONES. 

Como conclusion de todo lo anteriormente expuesto 
puede deducirse que: 

— El algoritmo estocastico de Kalman permite un ajuste 
riguroso del parametro a de la ecuacion que representa el 
proceso de filtracion a presion constante dando a la vez el 
grado de incertidumbre. 

— Es posible realizar el ajuste con un numero muy redu- 
cido de valores, tornados a lo largo de la prueba de laborato- 
rio lo cual, aunque no reduce la duration del experimento, 
ya que es necesario fmalizar la filtracion a fin de obtener el 
valor de la humedad final de la torta, permite un ahorro de 
tiempo cn la toma y tratamiento de los datos experimenla- 
les. 

NOMENCLATURA 

A : area de filtracion, cm 2 

a : coeficiente de la ecuacion f2) 

a : resistencia especifica de la torta, cm/g 

b : coeficiente de la ecuacion (2) 

k : constante de conversion de Newtons a Dinas 

e : eociente entre el tiempo y el volumcn filtrado, 

s/cm 3 


M : humedad final de la torta g. torta hinneda / g. torta 
seca. 

n : viscosidad del filtrado, g/cm.s. 

AP : presion total aplicada, Nw/cm 2 
Rm : resistencia del medio filtrante. cm’ 1 
p : densidad del filtrado g/cm 3 
t : tiempo de filtracion s 
V : volumen filtrado, cm 3 

x : concentration inicial de solidos, g solidos/g barbo- 
tina 

SV : covarianza del volumen 

ga : covarianza del parametro a 
ZV.a : covarianza cruzada de V y a. 
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